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ABSTRACT:  We  present  a  critique  of  Emanuel’s  steady-state  hurricane  model,  which  is  a  precursor  to  his  theory  for 
hurricane  potential  intensity  (PI).  We  show  that  a  major  deficiency  of  the  theory  is  the  tacit  assumption  of  gradient  wind 
balance  in  the  boundary  layer,  a  layer  that  owes  its  existence  to  gradient  wind  imbalance  in  the  radial  momentum  equation. 
If  a  more  complete  boundary-layer  formulation  is  included  using  the  gradient  wind  profiles  obtained  from  Emanuel’s 
theory,  the  tangential  wind  speed  in  the  boundary  layer  becomes  supergradient,  invalidating  the  assumption  of  gradient 
wind  balance.  We  show  that  the  degree  to  which  the  tangential  wind  is  supergradient  depends  on  the  assumed  boundary- 
layer  depth.  The  full  boundary-layer  solutions  require  a  knowledge  of  the  tangential  wind  profile  above  the  boundary  layer 
in  the  outer  region  where  there  is  subsidence  into  the  layer  and  they  depend  on  the  breadth  of  this  profile.  This  effect  is 
not  considered  in  Emanuel’s  theory.  We  argue  that  a  more  complete  theory  for  the  steady-state  hurricane  would  require 
the  radial  pressure  gradient  above  the  boundary  layer  to  be  prescribed  or  determined  independently  of  the  boundary  layer. 

The  issues  raised  herein  highlight  a  fundamental  problem  with  Emanuel’s  theory  for  PI,  since  that  theory  makes  the 
same  assumptions  as  in  the  steady-state  hurricane  model.  Our  current  findings  together  with  recent  studies  examining 
intense  hurricanes  suggest  a  way  forward  towards  a  more  consistent  theory  for  hurricane  PI.  Copyright  ©  2008  Royal 
Meteorological  Society 
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1.  Introduction 

In  the  first  of  what  has  turned  out  to  be  a  series 
of  very  influential  papers,  Emanuel  (1986,  henceforth 
E86)  presented  a  steady  axisymmetric  model  for  a 
mature  hurricane.  We  consider  this  paper  to  be  an 
important  milestone  in  tropical  cyclone  research  in  that 
it  re-focussed  attention  on  the  importance  of  the  radial 
gradient  of  sea  surface  moisture  fluxes  in  the  storm-scale 
energetics.  The  hurricane  model  described  therein  was  a 
prelude  to  the  development  of  an  axisymmetric  theory  for 
the  potential  intensity  (PI)  of  a  tropical  cyclone,  which 
we  refer  to  as  EPI-theory  (Emanuel,  1988;  Emanuel, 
1995;  Bister  and  Emanuel,  1998).  Since  its  inception, 
EPI-theory  has  been  called  upon  by  many  researchers  as 
a  standard  for  comparison  with  the  intensity  attained  in 
numerical  models  (e.g.  Frank  and  Ritchie,  2001;  Persing 
and  Montgomery,  2003)  or  assessments  of  possible 
changes  in  the  intensity  of  hurricanes  as  a  result  of  global 
warming  (e.g.  Knutson  and  Tuleya,  2004;  Emanuel,  2005; 
Bengtsson  et  al.,  2007).  At  the  present  time  it  appears  to 
be  the  only  such  theory  of  merit  for  these  applications 
(Camp  and  Montgomery,  2001).  Even  so,  there  are 
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indications  that  the  theory  is  deficient.  For  example, 
Persing  and  Montgomery  (2003)  have  shown  that  high- 
resolution  numerical  models  have  a  tendency  to  produce 
‘superintense’  storms,  superintense  meaning  that  they 
significantly  exceed  the  intensity  predicted  by  EPI-theory. 
Moreover,  the  calculated  PI  depends  sensitively  on  the 
assumed  relative  humidity  at  the  radius  of  maximum 
tangential  wind  speed,  which  Emanuel  generally  takes  to 
be  80%.  In  this  paper  we  draw  attention  to  a  fundamental 
inconsistency  of  the  hurricane  model  and  of  EPI-theory, 
namely  the  assumption  of  gradient  wind  balance  in 
the  boundary  layer,  both  inside  and  outside  the  radius 
of  maximum  tangential  wind  speed.  The  consequences 
of  this  assumption  for  Emanuel’s  hurricane  model  and 
EPI-theory  are  discussed  below  and  a  way  forward  is 
sketched. 

This  paper  is  organized  as  follows.  Section  2  provides 
a  brief  review  of  the  E86  hurricane  model,  with  the 
details  relegated  to  an  appendix.  Section  3  reviews  a 
more  complete  model  for  the  boundary  layer,  based  on 
the  work  of  Smith  (2003)  and  Smith  and  Vogl  (2008; 
henceforth  SV08),  and  examines  the  approximations 
made  by  Emanuel  in  terms  of  this  model.  Solutions  of 
the  more  complete  boundary-layer  model  with  Emanuel’s 
gradient  wind  profile  at  the  top  of  the  layer  are  presented 
and  discussed  in  Section  4.  Shown  also  in  this  section  is 
the  dependence  of  the  boundary-layer  flow  on  the  depth 
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of  the  layer  and  on  the  breadth  of  the  gradient  wind  profile 
at  its  top.  Section  5  discusses  some  implications  of  the 
calculations  and  Section  6  presents  the  conclusions. 


2.  The  E86  model  in  brief 

In  the  E86  model,  the  hurricane  vortex  is  assumed  to 
be  steady  and  circularly  symmetric  about  its  axis  of 
rotation.  The  boundary  layer  is  taken  to  have  uniform 
depth,  h,  and  is  divided  into  three  regions  as  shown 
in  Figure  1.  Regions  I  and  II  encompass  the  eye  and 
eyewall,  respectively,  while  Region  III  refers  to  that 
beyond  the  radius,  rm,  of  maximum  tangential  wind 
speed,  vm,  at  the  top  of  the  boundary  layer.  (Contrary 
to  Emanuel’s  assumption  in  this  figure,  observations 
show  that  rm  is  located  well  inside  the  outer  edge  of 
the  eyewall,  e.g.  Figure  3  of  Marks  et  al.  2008.  The 
significance  of  this  discrepancy  will  become  clearer  in 
section  5).  E86  takes  the  outer  radius  of  Region  II  to  be 
rm  on  the  basis  that  precipitation-driven  downdraughts 
may  be  important  outside  this  radius.  The  tangential 
wind  field  above  the  boundary  layer  is  assumed  to  be 
in  thermal  wind  balance  and  air  parcels  flowing  upwards 
and  outwards  into  the  upper  troposphere  are  assumed 
to  conserve  their  absolute  angular  momentum,  M ,  and 
saturation  moist  entropy,  s*  (calculated  reversibly).  These 
surfaces  are  assumed  to  flare  out  in  the  upper  troposphere. 
Here,  s*  is  defined  by: 

s*=cp\n9*,  (1) 


Figure  1.  Schematic  diagram  of  Emanuel's  1986  model  for  a  mature 
hurricane.  The  boundary  layer  is  assumed  to  have  constant  depth  h 
and  is  divided  into  three  regions  as  shown:  the  eye  (Region  I),  the 
eyewall  (Region  II)  and  outside  the  eyewall  (Region  III)  where  spiral 
rainbands  and  shallow  convection  emanate  into  the  vortex  above.  The 
absolute  angular  momentum  per  unit  mass,  M,  and  equivalent  potential 
temperature,  9e,  of  an  air  parcel  are  conserved  after  the  parcel  leaves 
the  boundary  layer  and  ascends  in  the  eyewall  cloud.  The  precise 
values  of  these  quantities  depend  on  the  radius  at  which  the  parcel  exits 
the  boundary  layer.  The  model  assumes  that  the  radius  of  maximum 
tangential  wind  speed,  rm,  is  located  at  the  outer  edge  of  the  eyewall 
cloud,  whereas  recent  observations  (e.g.  Figure  3  of  Marks  et  al.,  2008) 
indicate  that  it  is  closer  to  the  inner  edge.  This  figure  is  available  in 
colour  online  at  www.interscience.wiley.com/qj 


where  9*  is  the  reversible  saturation  equivalent  potential 
temperature  and  cp  denotes  the  specific  heat  at  constant 
pressure  of  dry  air.  Because  the  saturation  vapour  pres¬ 
sure  of  moist  air  is  a  unique  function  of  temperature  both 
s*  and  0*  are  state  variables.  E86  then  integrates  the 
thermal  wind  equation  upwards  along  these  surfaces  from 
radius  r  to  some  large  radius  rout  (>>  r)  to  obtain  a  rela¬ 
tionship  between  the  radial  rates  of  change  of  M  and  .v*  at 
the  top  of  the  boundary  layer,  z  =  h  (see  Equation  (A12) 
in  the  Appendix).  This  equation  may  be  further  integrated 
with  respect  to  radius  to  obtain  a  relationship  between  9* 
and  the  logarithm  of  the  Exner  function  at  the  top  of  the 
boundary  layer,  assuming  gradient  wind  balance  prevails 
at  this  height: 


where  TB  is  the  temperature  at  z  —  h,  Toul  is  the  tem¬ 
perature  on  the  v*  surface  at  rout  and  T  out  is  an  aver¬ 
age  of  this  temperature  weighted  with  the  saturation 
moist  entropy  of  the  outflow  angular  momentum  surfaces 
(Equation  (A  16)),  n  is  the  Exner  function,  /  is  the  Cori¬ 
olis  parameter  and  the  subscript  „  denotes  a  value  at  some 
large  radius  r  =  ra.  So  as  not  to  distract  from  the  main 
presentation,  the  details  are  given  in  the  appendix. 

The  flow  in  Regions  I  and  II  is  fully  determined 
by  a  simple  slab  formulation  for  the  boundary  layer 
from  which  a  second  functional  relationship  is  obtained 
between  M  and  s*  (Equation  (A17)).  The  two  relation¬ 
ships,  Equations  (A12)  and  (A17).  lead  inter  alia  to  an 
expression  for  the  tangential  wind  speed,  V,  at  z  —  h  in 
Region  II.  In  this  region  the  Rossby  number  is  large  com¬ 
pared  to  unity  and  the  Coriolis  term  may  be  neglected, 
giving 

liV2  =  ^CpiTv  -  Tout) (In  0*  -  In 9*),  at  z  =  h,  (3) 

where  9*s  is  the  saturation  equivalent  potential  tempera¬ 
ture  at  the  sea  surface  temperature,  Q  and  C/j  arc  sea 
surface  exchange  coefficients  for  enthalpy  and  momen¬ 
tum,  and  ji  —  VJV,  where  Vs  is  the  magnitude  of  the 
near-surface  wind.  Equation  (3)  states  that  in  Region  II, 
V  is  determined  locally  by  the  thermodynamic  disequi¬ 
librium  between  the  air  in  the  boundary  layer  and  the  sea 
surface  and  the  temperature  difference  between  the  top 
of  the  boundary  layer  and  the  outflow  temperature. 

E86’s  boundary-layer  formulation  in  Regions  I  and  II 
expresses  a  balance  between  radial  advection  and  surface 
gain  or  loss  of  azimuthal  momentum  and  specific  entropy. 
In  the  derivation  of  Equation  (3),  the  radial  velocity  is 
eliminated  so  that  the  formula  for  V 2  is  not  explicitly 
dependent  on  the  radial  component  of  velocity  in  the 
boundary  layer. 

Equations  (2)  and  (3)  lead  essentially  to  an  expres¬ 
sion  for  the  pressure  as  a  function  of  radius  (actually  the 


Copyright  ©  2008  Royal  Meteorological  Society 


Q.  J.  R.  Meteorol.  Soc.  134:  551-561  (2008) 
DOI:  10. 1002/qj 


HURRICANE  MODEL  AND  POTENTIAL  INTENSITY 


553 


logarithm  of  the  Exner  function)  at  the  top  of  the  bound¬ 
ary  layer  in  Regions  I  and  II  (Equations  (41)  and  (45)  of 
E86).  On  the  basis  that  precipitation-driven  downdraughts 
tend  to  offset  the  moistening  of  inflowing  air  parcels  in 
Region  III,  Emanuel  assumes  that  the  relative  humidity 
at  the  top  of  the  surface  layer  has  a  constant  value  of 
80%  all  the  way  inwards  to  rm,  an  assumption  that  is  not 
borne  out  by  observations  (e.g.  Figure  4d  of  Montgomery 
el  al.,  2006).  This  assumption  leads  to  a  second  equation 
relating  the  equivalent  potential  temperature  to  the  loga¬ 
rithm  of  the  Exner  function  and  the  relative  humidity  at 
the  top  of  the  surface  layer  (Equation  (A20)).  This  equa¬ 
tion,  when  combined  with  Equation  (2),  gives  an  expres¬ 
sion  for  the  logarithm  of  the  Exner  function  at  z  =  It 
in  Region  III  (Equation  (39)  of  E86).  With  the  assump¬ 
tion  of  gradient  wind  balance  at  z  —  h,  the  resulting  two 
equations  for  pressure  and  0*  at  this  level  completely 
determine  the  tangential  wind  speed  at  the  top  of  the 
boundary  layer  at  all  radii. 

Note  that  the  tangential  wind  speed  at  the  top  of 
Region  III  is  obtained  only  from  thermodynamic  con¬ 
siderations  in  the  boundary  layer;  the  dynamics  of  the 
boundary  layer  are  completely  ignored.  It  will  be  argued 
below  that  the  tacit  assumption  of  gradient  wind  balance 
in  the  boundary  layer  in  Regions  I  and  II  and  the  neglect 
of  boundary-layer  dynamics  in  Region  III  represent  a  fun¬ 
damental  limitation  of  Emanuel’s  theory  and  leads  to  an 
inconsistency  with  important  ramifications. 


3.  The  boundary  layer 

To  put  E86’s  assumptions  regarding  the  boundary  layer 
in  perspective,  we  review  first  the  simple,  but  more  com¬ 
plete  model  of  the  boundary  layer  of  a  steady  axisym- 
rnetric  hurricane-like  vortex  on  an  /-plane  developed  by 
Smith  (2003)  and  SV08.  We  examine  then  the  conse¬ 
quences  of  assuming  gradient  wind  balance  in  the  layer. 

3.1.  A  slab  model  for  the  boundary  layer 

The  boundary  layer  in  SV08  is  assumed  again  to  have 
uniform  depth,  h,  and  constant  density.  (SV08  considered 
also  the  variable  depth  case,  but  for  simplicity  the  focus 
here  is  on  the  constant-depth  boundary  layer  assumed 
by  E86.)  In  a  cylindrical  coordinate  system  (r,<f>,z), 
the  vertically  integrated  equations  expressing  the  local 
budgets  of  radial  momentum,  azimuthal  momentum,  heat 
or  moisture,  and  mass  continuity  can  be  written  in  the 
following  form: 


dwb  Wh_  +  W sc  (Ugr  “  %) 

Ufa—  —  -  U\) 

dr  h  r 

-  /Ogr  -  %)  -  —J—  (Mb  + 

h 

dub  u’i,  +  wsc  .Vb 

Mb—  =  - ; - Ob  -  fgr)  -  (—  +  /)wb 
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Cd. 
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- r^Ob  +  vl)lvb, 
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r  h 


(6) 

(7) 


where  Mb  and  nb  are  the  radial  and  azimuthal  components 
of  velocity  in  the  boundary  layer,  ugr(r)  and  ?/,’/,  are 
the  tangential  velocity  and  vertical  velocity  at  the  top 
of  the  boundary  layer,  Wh- =  \(wh  —  |t%|),  Xb  is  a 
scalar  quantity,  which  could  be  the  dry  static  energy, 
the  specific  humidity,  or  the  specific  entropy,  /  is  the 
Coriolis  parameter,  Co  is  the  surface  drag  coefficient, 
Cx  is  the  surface  transfer  coefficient  for  Xb>  X/;+  is 
the  value  of  /  just  above  the  boundary  layer,  and 
Xs  is  the  value  of  /  at  the  sea  surface.  The  terms 
involving  wsc  represent  turbulent  fluxes  at  the  top  of 
the  boundary  layer  (arising  from  rainbands,  shallow 
convection,  or  smaller-scale  turbulent  structures)  and 
the  term  Xb  represents  the  effects  of  radiative  cooling 
and  dissipative  heating  when  Xb  is  taken  to  be  the  dry 
static  energy  or  specific  entropy.  Consistent  with  the  slab 
boundary-layer  formulation,  the  quantities  ub,  ub  and  /b 
are  assumed  to  be  independent  of  depth.  Note  that  w;,_  is 
nonzero  only  when  w/,  <  0,  in  which  case  it  is  equal  to 
wh.  Thus  the  terms  involving  w/,_  represent  the  transport 
of  properties  from  above  the  boundary  layer  that  may  be 
different  from  those  inside  the  boundary  layer.  For  the 
calculations  presented  in  sections  4.1  and  4.2  we  take 
Co  to  be  a  constant,  equal  to  2.0  x  10-3,  the  value  use 
by  E86.  For  those  in  section  4.3,  we  follow  SV08  and 
take  Co  =  Coo  +  Coi|ub|,  where  Coo  =  0.7  x  10~3  and 
Coi  =  6.5  x  10-5  for  wind  speeds  less  than  20  m  s_1 
and  Co  =  2.0  x  10~3,  a  constant,  for  larger  wind  speeds. 
These  values  are  based  on  our  interpretation  of  Figure  5 
from  Black  et  al.  (2007).  In  the  calculations  described  in 
section  4,  we  consider  only  dynamical  effects,  so  that  a 
value  for  Cx  is  not  required. 

Substitution  of  Equation  (7)  into  (4)  gives  an  expres¬ 
sion  for  in  i, : 


wh 


1  +  a 


1 

Mb 


gr 


+  f(V gr  -  fb) 


+  -f~(Ub  +  Mb)  2  Mb  )  -  — 
n  r 


(8) 


where  a  is  zero  if  the  expression  in  square  brackets  is 
positive  and  unity  if  it  is  negative.  With  this  expression 
for  Wh,  Equations  (4)— (6)  together  with  (8)  form  a  system 
of  ordinary  differential  equations  that  may  be  integrated 
radially  inwards  from  some  large  radius  R  to  find  «b,  nb 
and  Xb  as  functions  of  r,  given  values  of  these  quantities 
at  r  =  R  as  well  as  the  radial  profile  vp  (r). 


3.2.  E86’s  approximations  for  the  boundary  layer 

Emanuel  writes  Equation  (5)  in  terms  of  the  absolute 
angular  momentum  in  the  boundary  layer,  Mb  —  n>b  + 
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fr2 /2,  and  approximates  this  equation  in  Region  II, 
where  Wh  >  0,  as 


in  Equation  (5)  can  be  reasonably  approximated  by  vgI, 
the  equation  predicts  that 


U  b 


d  Mb 

dr 


(9) 


Here  it  is  assumed  that  wsc  =  0  and  that  Mb  <<  ty  in 
the  drag  term.  Note  that,  in  general,  knowledge  of  Mb  is 
required  for  the  determination  of  Mb.  However,  Emanuel 
does  not  use  the  radial  momentum  equation  to  determine 
«b,  as  his  main  focus  is  to  obtain  an  expression  relating 
the  specific  entropy,  sb,  to  Mb  (Equation  (A17)).  (In  fact, 
E86  uses  Equation  (9)  to  determine  Mb,  having  obtained 
the  radial  pressure  distribution  through  his  Equations  (39) 
and  (41)  and  having  tacitly  assumed  gradient  wind 
balance  to  obtain  t>b.)  In  the  region  where  w/,  >  0,  the 
equation  for  the  specific  entropy  is 


Ck 

h 


VbCs*  -^b), 


(10) 


where  s*  is  the  saturation  specific  humidity  at  the  sea 
surface  and  again  it  is  assumed  that  wsc  —  0  and  the  total 
wind  speed  has  been  approximated  by  the  tangential  wind 
speed. 

E86’s  assumption  that  air  leaving  the  boundary  layer 
conserves  its  absolute  angular  momentum  implies  that 
ugr  —  t’l,  where  Wh  >  0.  The  assumption  also  that  ugr  is 
in  gradient  wind  balance  implies  that  vb  is  in  gradient 
wind  balance.  This  is  a  rather  strong  assumption  for  the 
boundary  layer  in  the  inner  core  of  a  rapidly  rotating 
vortex  and  although  it  has  been  made  by  previous  authors 
(e.g.  Ooyarna,  1969),  we  are  not  aware  of  any  rigorous 
justification  for  it.  In  fact  it  is  not  supported  by  a  scale 
analysis  of  the  boundary-layer  equations  (e.g.  Smith, 
1968).  Ooyarna  was  certainly  aware  of  the  limitations  of 
the  assumption  and  wrote  in  an  unpublished  manuscript 
in  1968: 


...  it  appears  that  the  weakest  hypothesis  in  [ his ]  original 
model  is  the  use  of  the  balance  approximation  in  the 
boundary  layer. 

In  this  manuscript,  Ooyarna  went  on  to  show  that 
the  solutions  in  a  calculation  with  a  more  complete 
boundary-layer  formulation  were  more  realistic  than 
those  with  a  balanced  boundary-layer  formulation.  Indeed 
it  is  precisely  the  lack  of  gradient  wind  balance  in  the 
boundary  layer  that  gives  rise  to  the  ‘frictionally  driven' 
inflow  in  the  layer. 

While  inflow  is  theoretically  possible  in  a  boundary 
layer  that  is  in  approximate  gradient  wind  balance, 
the  balance  assumption  can  be  justified  only  if  the 
radial  acceleration  and  radial  friction  terms  are  small 
compared  with  the  radial  pressure  gradient  and  the 
sum  of  the  centrifugal  and  Coriolis  forces.  In  such  a 
‘balanced’  formulation,  the  radial  flow  is  determined  by 
the  (sic)  tangential  momentum  equation.  With  Emanuel’s 
assumption  that  the  total  wind  speed  in  the  friction  term 


Ub  —  -CV  gr,  (11) 

where  c  —  Qugr/(It£a),  and  fa(=  £  +  /)  and  f  are  the 
absolute  vorticity  and  relative  vorticity  of  the  gradient 
wind,  t'g,.,  respectively.  Other  processes  could  contribute 
also  to  radial  motion  in  a  boundary  layer  that  is  closely  in 
gradient  wind  balance.  One  example  is  a  radial  buoyancy 
gradient  above  the  boundary  layer  associated  with  moist 
convective  processes  (e.g.  Smith,  2000;  Smith  et  al., 
2005). 

In  the  next  section  we  examine  solutions  of  the  dynam¬ 
ical  component  of  the  full  boundary-layer  equations  (4), 
(5)  and  (7)  with  the  gradient  wind  speed  ugr  obtained  by 
E86.  We  show  that  these  solutions  are  incompatible  with 
the  assumption  in  the  E86  model  that  ugr  =  nb  where 
Wh  >  0.  We  show  further  that  the  lack  of  any  dynamical 
constraint  in  the  boundary  layer  in  Region  III  other  than 
the  tacit  assumption  of  gradient  wind  balance  is  another 
major  deficiency  of  the  theory. 


4.  Calculations 

4.1.  The  E8 6  gradient  wind  profile 

Figure  2  shows  calculations  of  the  full  boundary-layer 
equations  of  Section  3.1,  taking  the  gradient  wind  speed 
profile  vgr(r)  and  other  parameters  the  same  as  those 
obtained  by  E86.  In  particular  /  =  6.83  x  10~5  s-1, 
corresponding  with  a  latitude  of  28  °N,  h  —  1000  m, 
CD  =  2.0  x  1(T3,  Ts  =  27  °C,  TB  =  27  °C  and  Tout  = 
—67  °C.  The  radial  profile  of  ugr  is  obtained  by  solving  the 
gradient  wind  equation  with  the  pressure  profile  derived 
from  the  pair  of  expressions  relating  lnjr  and  9t*  in 
E86,  namely  Equations  (39)  and  (41),  using  the  parameter 
values  detailed  in  that  paper.  The  integration  in  the  full 
boundary-layer  calculation  starts  at  a  radius  of  375  km, 
where  the  gradient  wind  speed  (only  1.73  m  s_1)  is  small 
enough  to  justify  the  neglect  of  the  nonlinear  acceleration 
terms  in  the  equations  (section  4  of  Smith,  2003).  (Note 
that  beyond  a  radius  of  400  km,  the  tangential  wind 
in  Emanuel’s  calculation  is  anticyclonic  and  just  inside 
this  radius,  at  about  396  km,  the  profile  is  inertially 
unstable.)  Figure  2(a)  compares  the  full  solution  for  the 
tangential  wind  speed  in  the  boundary  layer,  ub,  with  the 
imposed  gradient  wind  speed  ugr.  It  compares  also  the  full 
solution  for  the  radial  wind  speed,  Mb,  with  that  obtained 
from  Equation  (11)  based  on  the  balance  assumption  that 
ngr  =  Vb  as  made  by  E86,  and  assuming  that  wsc  =  0. 
We  designate  the  balanced  solution  for  Mb  as  uE  and 
that  for  the  corresponding  vertical  motion  at  the  top  of 
the  boundary  layer  as  wE.  We  calculate  the  latter  from 
the  continuity  equation  (7)  using  centred  differences.  The 
profiles  of  vertical  velocity  at  the  top  of  the  boundary 
layer  in  the  full  solution,  u>i, ,  is  compared  with  that  in 
the  balanced  solution  in  Figure  2(b).  It  is  worth  noting  at 
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Figure  2.  (a)  Radial  profiles  of  boundary-layer  radial  (Mb)  and  tangential  (1%)  wind  components  and  the  total  wind  speed  y/ («b2  +  rb2)  (denoted 
vv)  from  the  full  boundary-layer  solution,  and  the  tangential  wind  speed  above  the  boundary  layer  (Vgr)  as  obtained  by  E86  (solid  curve).  Also 
shown  is  the  radial  flow  obtained  from  the  balanced  solution  determined  from  Equation  (11)  and  denoted  by  ue .  Units  are  m  s_1.  For  plotting 
convenience,  the  signs  of  Mb  and  ue  have  been  reversed.  The  profile  of  vgr  is  indicated  by  the  unmarked  solid  curve,  (b)  Corresponding  radial 
profiles  of  vertical  velocity  at  the  top  of  the  boundary  layer  (ro/,)  and  that  in  the  balanced  solution  (u>e)-  Units  are  cm  s-1 .  The  thin  vertical  line 
in  (a)  and  (b)  marks  the  radius  of  maximum  vgr,  the  boundary  between  Regions  II  and  III  in  Figure  1.  (c)  Radial  profile  of  the  coefficient  c  in 
Equation  (11).  (d)  Radial  profiles  of  the  three  terms  on  the  right-hand  side  of  the  radial  momentum  equation,  Equation  (4),  and  their  sum  (the 
solid  line)  for  the  full  solution.  The  designation  ‘wu\  ‘nif  and  ‘fri’  refer  to  the  first,  second  and  third  terms  in  the  equation,  representing  the 
downward  advection  of  radial  momentum,  the  net  inward  force  and  the  frictional  force,  respectively.  This  figure  is  available  in  colour  online  at 

www  .interscience .  wiley  .com/qj 


this  point  that  this  balanced  solution  agrees  closely  with 
that  shown  by  E86  in  his  Figure  12. 

In  the  full  and  balanced  calculations,  the  radial  wind 
component  increases  inwards  to  a  certain  radius  and  then 
decreases.  However,  there  are  significant  quantitative  dif¬ 
ferences  in  the  profiles.  In  the  balanced  solution,  the 
maximum  inflow  of  about  12  m  s~ 1  occurs  at  a  compar¬ 
atively  large  radius  (130  km),  while  in  the  full  solution  it 
occurs  at  52  km,  a  little  outside  the  radius  of  maximum 
gradient  wind  speed  (35.8  km).  These  differences  occur 
despite  the  fact  that  beyond  100  km  in  radius,  ub  is  at 
most  18%  smaller  than  ngr,  showing  that  the  degree  of 
gradient  wind  imbalance  is  important.  The  decline  in  ue 
from  such  a  large  radius  is  a  result  of  the  decline  in  the 


parameter  c  with  decreasing  radius  (Figure  2(c)),  which 
is  larger  than  the  rate  at  which  ngr  increases.  The  discon¬ 
tinuity  in  ue  at  r  =  rm  is  a  result  of  the  discontinuity  of 
the  relative  vorticity  £  at  this  radius,  which  leads  to  a 
discontinuity  in  c.  As  expected  there  are  correspondingly 
large  differences  in  the  profiles  of  vertical  velocity  at  the 
top  of  the  boundary  layer  (Figure  2(b)).  In  particular,  the 
change  from  descent  at  large  radii  to  ascent  at  small  radii 
occurs  at  a  much  smaller  radius  in  the  full  calculation: 
107  km  compared  with  230  km. 

Of  particular  significance  is  the  difference  between  ub 
and  ngr  in  the  inner  core  region,  near  the  radius  of  max¬ 
imum  gradient  wind  speed.  Here  the  tangential  wind  in 
the  boundary  layer  becomes  supergradient  (i.e.  i>b  exceeds 


Copyright  ©  2008  Royal  Meteorological  Society 


Q.  J.  R.  Meteorol.  Soc.  134:  551-561  (2008) 
DOI:  10. 1002/qj 


556 


R.  K.  SMITH  ETAL., 


Dgr),  which  is  incompatible  with  Emanuel’s  assumption 
that  vgr  is  equal  to  i>b  at  radii  where  w/,  >  0.  In  other 
words,  Emanuel’s  calculated  PI  (i.e.  vm)  is  exceeded 
when  a  more  complete  boundary  layer  formulation  is 
used.  The  occurrence  of  supergradient  winds  is  a  reflec¬ 
tion  of  the  strong  radial  inflow  which  advects  absolute 
angular  momentum  at  a  rate  larger  than  it  can  be  removed 
locally  by  the  frictional  torque  (SV08).  As  soon  as  the 
tangential  wind  speed  becomes  supergradient,  all  forces 
in  the  radial  momentum  equation  act  outwards  and  lead 
to  a  rapid  deceleration  of  the  inflow.  In  the  full  boundary- 
layer  solution,  the  radial  flow  becomes  zero  at  some  finite 
radius  and  the  boundary-layer  model  becomes  singular  at 
this  radius.  In  reality  we  would  expect  the  inflow  to  be 
expelled  upwards  before  this  radius,  carrying  its  hori¬ 
zontal  momentum  with  it.  If  the  upflow  remains  out  of 
balance,  we  would  expect  it  to  flow  outwards  imme¬ 
diately  above  the  inflow  layer,  a  behaviour  which  is 
shown  by  full  numerical  solutions  (e.g.  Figures  3c,  6c 
of  Montgomery  et  al.,  2001;  Figure  lb  of  Persing  and 
Montgomery,  2003). 

Figure  2(d)  shows  the  radial  variation  of  the  force 
terms  in  the  radial  momentum  equation.  Equation  (4). 
The  term  representing  the  downward  transport  of  radial 
momentum,  that  proportional  to  Wh~,  is  non-zero  only 
in  the  outer  region  and  is  small  compared  with  the 
other  terms.  At  larger  radii,  the  net  inward  force  (the 
difference  between  the  inward  pressure  gradient  and 
outward  centrifugal  and  Coriolis  forces)  is  larger  in 
magnitude  than  the  outward  frictional  force.  Moreover, 
the  inward  radial  acceleration,  which  is  equal  to  the  net 
total  radially  inward  force,  is  particularly  large  at  radii 
less  than  150  km. 

4.2.  Dependence  on  boundary-layer  depth 

The  boundary-layer  constraint  in  Emanuel’s  theory  is 
independent  of  the  assumed  boundary-layer  depth.  (Note 
that  the  depth  cancels  in  applying  E86’s  boundary-layer 


Depth  —  600  m 


Figure  3.  As  Figure  2(a).  but  for  a  boundary  layer  depth  of  600  m.  This 
figure  is  available  in  colour  online  at  www.interscience.wiley.com/qj 


formulation  to  derive  Equation  (A17)  in  the  appendix.) 
However,  this  depth  has  a  significant  influence  on  the 
full  boundary-layer  solution  because  the  effective  drag 
in  the  boundary  layer  is  inversely  proportional  to  the 
depth  (SV08).  For  this  reason  we  repeated  the  foregoing 
boundary-layer  calculations  for  a  boundary-layer  depth 
of  600  m.  These  calculations  are  shown  in  Figure  3.  The 
increased  effective  friction  leads  to  a  larger  reduction  of 
the  tangential  wind  speed  in  the  boundary  layer  than 
in  the  earlier  calculation  and  therefore  to  a  larger  net 
inward  force  and  a  larger  inward  acceleration.  Conse¬ 
quently  the  maximum  inflow  is  considerably  larger  than 
before  (36  m  s~'  instead  of  19  ms"1)  and  occurs  at  a 
smaller  radius  (32  km  instead  of  52  km).  On  the  other 
hand,  the  balanced  solution  changes  only  in  magnitude 
and  not  in  shape,  whereupon  the  maximum  occurs  at 
130  km  as  before.  This  result  follows  directly  from  Equa¬ 
tion  (11)  because  the  decreased  depth  simply  increases 
the  coefficient  c  by  a  constant  factor  at  all  radii  and  the 
gradient  wind  profile  is  the  same.  The  fact  that  the  maxi¬ 
mum  tangential  wind  speed  in  the  boundary  layer  in  this 
calculation  is  considerably  higher  than  in  the  previous 
one  implies  that  the  PI  of  the  steady  vortex  is  sensi¬ 
tive  to  the  boundary-layer  depth,  an  important  point  not 
emphasized  in  E86  and  his  subsequent  papers.  (Whereas 
the  E86  model  and  the  more  complete  boundary-layer 
model  furnish  non-negligible,  but  modest  differences  in 
the  maximum  tangential  wind  (~10-20%),  it  should  be 
remembered  that  the  boundary-layer  model  used  here 
precludes  any  thermodynamic  and  dynamic  feedbacks 
between  the  boundary  layer  and  interior  flow.  For  sev¬ 
eral  reasons,  this  feedback  is  thought  to  be  quantitatively 
significant;  see  section  5). 

4.3.  Dependence  on  vortex  size 

Given  the  importance  of  the  radial  acceleration  in  the 
boundary  layer  as  demonstrated  above,  the  inclusion  of 
boundary-layer  dynamics  in  Region  III  of  Emanuel’s 
model  may  be  expected  to  have  important  consequences 
for  the  tangential  wind  maximum  also.  We  illustrate  these 
consequences  by  a  third  set  of  calculations  to  emphasize 
the  dependence  of  the  maximum  boundary-layer  wind 
speed  on  the  vortex  size.  These  calculations  are  based 
on  solutions  of  the  full  boundary-layer  equations  with 
the  different  profiles  of  gradient  wind  speed  shown  in 
Figure  4.  These  profiles  are  defined  in  Smith  (2003) 
and  are  inertially  stable  for  the  value  of  /  used  earlier. 
The  solutions  for  these  profiles  are  shown  in  Figure  5 
for  a  boundary-layer  depth  of  800  m,  a  radially  varying 
drag  coefficient  Co  as  discussed  in  Section  3.1,  and 
with  u;Sc  =  —5.7  cm  s-1,  the  value  used  in  SV08.  Note 
that  there  is  a  clear  dependence  of  the  solution  on 
storm  size,  as  might  be  characterized,  for  example,  by 
the  radius  of  gale-force  winds  (17  m  s_1)  above  the 
boundary  layer.  As  the  storm  size  decreases,  the  radius 
of  maximum  inflow  decreases  and  the  maximum  inflow 
increases.  Moreover,  the  radius  at  which  the  vertical 
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Figure  4.  Four  radial  profiles  of  tangential  wind  speed,  ugr(r),  at 
the  top  of  the  boundary  layer  used  for  the  calculations  shown 
in  Figure  5.  The  thin  horizontal  line  indicates  the  radius  of  gale 
force  winds  (17  m  s-1).  This  figure  is  available  in  colour  online  at 
www.interscience.wiley.com/qj 

velocity  changes  sign  decreases  (figure  not  shown).  To 
the  extent  that  the  intensity  is  controlled  by  boundary- 
layer  dynamics,  these  solutions  show  a  clear  dependence 
on  the  size  of  the  outer  circulation,  so  that  the  PI  of 
midget  storms  may  be  expected  to  be  different  from 
that  of  broad  storms.  These  solutions  highlight  the 
dependence  of  the  flow  at  all  radii  in  the  boundary  layer 
on  the  size  of  the  vortex  above. 

5.  Discussion 

Using  the  gradient  wind  profile  predicted  by  Emanuel’s 
steady-state  hurricane  model  in  conjunction  with  a  more 
complete  formulation  of  the  boundary  layer  generally 


leads  to  the  occurrence  of  supergradient  winds  in  the 
boundary  layer  in  the  high  wind  region  of  the  vor¬ 
tex.  These  are  incompatible  with  a  key  assumption  in 
Emanuel’s  derivation  of  the  gradient  wind  profile  that 
requires  it  to  be  equal  to  that  in  the  boundary  layer 
where  the  flow  is  upwards  out  of  the  boundary  layer. 
Moreover,  the  degree  to  which  the  boundary-layer  winds 
are  supergradient  increases  as  the  boundary-layer  depth 
decreases.  In  reality,  the  vertical  advection  of  the  super¬ 
gradient  winds  out  of  the  boundary  layer  would  lead  to 
outflow  until  a  radius  is  achieved  at  which  the  pressure 
gradient  is  matched  to  that  which  can  be  sustained  by 
the  mass  distribution.  Of  course,  this  effect  cannot  be 
captured  by  a  one-layer  model,  but  it  is  significant  that 
calculations  in  which  the  boundary  layer  is  allowed  to 
adjust  to  an  outer  flow  do  show  such  behaviour  (e.g. 
Figures  3c,  6c  of  Montgomery  et  al.,  2001;  Figure  lb  of 
Persing  and  Montgomery,  2003). 

The  dependence  of  the  radius  at  which  subsidence 
at  large  radii  changes  to  ascent,  rup,  as  well  as  the 
predicted  radial  profiles  of  «b,  ub  and  Wh,  on  the 
tangential  wind  profile  above  the  boundary  layer  where 
there  is  subsidence  into  it  shows  that  the  dynamics  of 
the  boundary  layer  in  Region  HI  of  Figure  1  cannot  be 
ignored. 

The  foregoing  considerations  suggest  an  alternative 
subdivision  of  the  boundary  layer  to  that  in  Figure  1. 
This  alternative  is  sketched  in  Figure  6  and  is  based 
on  whether  the  top  of  the  boundary  layer  is  an  inflow 
boundary  (Region  B,  r  >  rup)  or  an  outflow  boundary 
(Region  A,  r  <  rup).  In  Region  B,  the  boundary  layer  is 
directly  influenced  by  the  vortex  above  through  the  radial 
pressure  gradient  at  the  top  of  the  layer  and  through  the 
downward  advection  of  free  vortex  properties  such  as 
moisture,  heat  and  momentum.  Except  possibly  through 
the  occurrence  of  moist  convection,  there  is  no  essen¬ 
tial  feedback  to  the  free  vortex.  (An  important  excep¬ 
tion  arises  with  the  occurrence  of  spiral  rainbands  and 


Tangential  wind  profiles 


delta  800  m  delta  800  m 


Figure  5.  Radial  profiles  of  (a)  the  radially  inward  and  (b)  the  tangential  components  of  wind  speed  in  the  boundary  layer  for  the  four  vortex 
profiles  shown  in  Figure  4.  This  figure  is  available  in  colour  online  at  www.interscience.wiley.com/qj 
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Figure  6.  Modified  conceptual  model  of  the  hurricane  inner-core  region 
motivated  by  the  findings  herein  together  with  recent  observational  and 
modelling  studies.  Air  subsides  into  the  boundary  layer  for  r  >  ru p 
and  ascends  out  of  the  boundary  layer  for  r  <  rUp-  The  frictionally 
induced  net  inward  force  in  the  outer  region  produces  a  radially 
inward  jet  at  r  =  rup.  The  subsequent  evolution  of  this  jet  depends 
on  the  bulk  radial  pressure  gradient  that  can  be  sustained  by  the 
mass  distribution  at  the  top  of  the  boundary  layer.  The  jet  eventually 
generates  supergradient  tangential  winds  whereafter  the  radial  flow 
rapidly  decelerates  and  turns  upwards  and  outwards.  When  the  outflow 
has  adjusted  to  the  radial  pressure  gradient  that  is  sustained  by  the  mass 
field,  the  flow  turns  upwards  into  the  eyewall  clouds.  See  sections  5 
and  6  for  further  details.  This  figure  is  available  in  colour  online  at 
www.interscience.wiley.com/qj 

the  corresponding  formation  of  one  or  more  secondary 
eyewalls  (Houze  et  al.,  2007;  Terwey  and  Montgomery, 
2008).  These  asymmetric  processes,  their  coupling  to  the 
boundary  layer  and  the  free  axisymmetric  vortex  are 
not  yet  well  understood  and  consequently  lie  beyond 
the  scope  of  the  present  model.)  However,  in  Region  A, 
boundary-layer  properties  are  advected  into  the  free  vor¬ 
tex  and  have  a  profound  influence  on  its  structure.  We 
may  think  of  the  boundary-layer  flow  in  Region  B  as 
producing  an  inward  radial  jet  at  r  —  rup,  the  strength  of 
which  depends  on  the  gradient  wind  profile  at  larger  radii 
as  well  as  on  the  boundary-layer  depth.  The  boundary- 
layer  dynamics  in  Region  A  determine  the  fate  of  this  jet, 
but  the  details  depend  inter  alia  on  the  radial  pressure 
gradient  at  the  top  of  the  boundary  layer,  i.e.  there  is  a 
substantial  two-way  feedback  between  the  boundary  layer 
and  the  free  vortex  in  this  region.  These  details  depend 
also  on  the  boundary-layer  depth.  The  radial  pressure  gra¬ 
dient  in  the  boundary  layer  is  probably  still  determined  in 
large  measure  by  the  mass  distribution  in  the  free  vortex, 
with  possible  exceptions  in  localized  regions  near  where 
inflow  turns  to  upflow  and  possibly  outflow  (see  below). 
However  the  free  vortex  can  be  expected  to  be  strongly 
influenced  by  the  radial  distribution  of  mass,  momentum 
and  moisture  that  leave  the  boundary  layer. 

The  calculations  presented  here,  supported  by  those  of 
SV08,  show  that  the  tangential  winds  tend  to  become 
supergradient  in  the  inner  core  and,  as  a  result,  the  radial 
flow  rapidly  decelerates  until  the  tangential  component 
becomes  subgradient  again,  or  the  radial  wind  becomes 
zero  (a  point  at  which  the  boundary-layer  equations 
in  a  one-layer  model  become  singular  and  a  more 


sophisticated  technique  beyond  the  scope  of  this  study 
is  required  for  matching  the  solutions  inside  and  outside 
this  radius).  In  either  case,  the  flow  out  of  the  boundary 
layer  increases  markedly.  If  the  winds  carried  upwards 
retain  their  supergradient  character  they  will  surely  flow 
with  a  significant  component  outwards  until  they  have 
come  into  gradient  wind  adjustment  with  the  mass  field 
aloft.  At  this  point  they  would  be  expected  to  turn 
upwards  into  the  eyewall.  While  parts  of  this  scenario 
are  speculative  at  this  stage,  the  foregoing  ideas  would 
explain  the  observations  of  a  skirt  of  moderate  to  high 
radar  reflectivity  adjacent  to  the  main  eyewall  (e.g. 
Figures  5-7  of  Aberson  et  al.,  2006;  Figure  3  of  Marks 
et  al.,  2008),  but  still  within  the  ‘visible’  eye  defined  by 
the  upper-tropospheric  boundary  of  clear  and  cloudy  air 
seen  in  high-resolution  satellite  imagery  (e.g.  Figure  2  of 
Bell  and  Montgomery,  2008)  and  they  are  consistent  with 
the  calculations  of  Montgomery  et  al.  (2001,  Figures  3c, 
6c)  and  Persing  and  Montgomery  (2003,  Figure  lb). 

Within  the  context  of  the  axisymmetric  model,  the  ther¬ 
modynamic  consequences  of  the  overshoot/adjustment 
region  have  been  demonstrated  to  be  non-trivial  as  moist 
air  near  the  surface  and  within  the  maximum  tangential 
wind  (including  the  outer  part  of  the  ‘eye’)  can  be  drawn 
into  the  main  eyewall  above  the  shallow  inflow  layer. 
This  low-level  air  generally  possesses  higher  equivalent 
potential  temperature  than  air  found  at  the  radius  of  max¬ 
imum  wind  due  to  a  lower  surface  pressure  and  non-zero 
surface  winds  and  contributes  additional  heat  and  local 
buoyancy  to  the  eyewall  (Persing  and  Montgomery,  2003; 
Cram  et  al.,  2007).  The  net  result  is  an  enhancement 
of  the  radial  gradient  of  equivalent  potential  tempera¬ 
ture  above  the  inflow  layer  that  supports  strong  tangential 
winds  in  accordance  with  axisymmetric  thermal  wind  bal¬ 
ance  above  the  boundary  layer  (Appendix  in  Montgomery 
et  al.,  2006).  In  light  of  these  findings,  together  with  the 
recognition  that  shear  instability  and  coherent  vortex  sub¬ 
structures  bordering  the  eye  and  eyewall  will  contribute  to 
the  aforementioned  adjustment  process  (Schubert  et  al., 
1999;  Montgomery  et  al.,  2002;  Braun  et  al.,  2006),  we 
believe  that  both  the  initial  vortex  structure  and  interac¬ 
tions  between  the  eye  and  eyewall  region  are  important 
elements  of  intense  storms  and  should  be  accounted  for 
in  hurricane  intensity  theory. 

Note  that  much  of  the  foregoing  discussion  was  based 
on  the  assumption  that  the  main  dynamical  processes 
in  tropical  cyclones  are  axisymmetric.  However,  recent 
calculations  by  Nguyen  et  al.  (2008)  have  highlighted  the 
importance  of  asymmetric  processes  in  the  intensification 
of  these  storms.  A  comprehensive  synthesis  of  these 
findings  with  the  insights  obtained  here  is  a  goal  for  future 
work. 


6.  Conclusions 

We  have  shown  that  the  tacit  assumption  of  gradient 
wind  balance  in  the  boundary  layer  is  a  major  deficiency 
of  Emanuel’s  steady-state  hurricane  model  and  also. 
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by  implication,  his  theory  for  the  potential  intensity  of 
hurricanes.  Although  the  vertically  integrated  tangential 
wind  in  the  boundary  layer  is  usually  no  more  than 
15-20%  less  than  its  gradient  wind  counterpart,  a  fact 
that  makes  gradient  wind  balance  a  seemingly  defensible 
zero-order  approximation  locally,  we  have  shown  that  the 
global  consequences  of  this  simplification  on  the  inner- 
core  structure  of  intense  storms  are  non-trivial.  Indeed, 
the  boundary  layer  owes  its  existence  to  gradient  wind 
imbalance  that  results  from  a  reduction  of  the  tangential 
wind  speed  by  friction.  When  such  imbalance  is  allowed 
for  by  the  inclusion  of  a  non-trivial  radial  momentum 
equation  in  the  theory,  the  boundary-layer  flow  depends 
on  the  tangential  wind  structure  above  the  boundary  layer, 
a  feature  that  must  be  taken  into  account  in  an  improved 
theory  for  hurricane  PI. 

We  conclude  that  it  is  not  permissible  to  make  the 
gradient  balance  assumption  in  the  inner  region  and 
that,  in  a  realistic  model  of  a  hurricane,  the  radial 
pressure  gradient  above  the  boundary  layer  must  be 
prescribed  or  determined  independently  of  the  boundary 
layer.  Nevertheless,  even  in  this  case,  the  solutions  show 
a  mismatch  between  the  predicted  mean  winds  in  the 
boundary  layer  and  those  prescribed  above  where  the 
flow  is  out  of  the  layer.  This  mismatch  suggests  that 
the  outflow  jet  found  above  the  inflow  layer  in  full 
numerical  solutions  for  the  boundary  layer,  together  with 
the  flow  above  it,  is  a  means  by  which  the  flow  exiting 
the  boundary  layer  adjusts  to  the  radial  pressure  gradient 
associated  with  the  vortex  above  the  boundary  layer.  The 
implication  would  be  that  a  more  complete  formulation 
of  the  (steady)  boundary  layer  in  the  inner  core  region  of 
a  tropical  cyclone  using  a  slab-type  formulation  would 
require  at  least  two  layers  including  one  to  represent  the 
outflow  jet.  This  layer  is  required  to  allow  the  radial  and 
tangential  wind  fields  to  adjust  to  the  radial  pressure 
gradient  implied  by  the  mass  distribution  in  the  free 
troposphere.  Such  a  formulation  would  appear  to  be  a 
necessary  component  of  a  more  consistent  and  accurate 
theory  for  hurricane  PI  and  such  a  theory  must  take  into 
account  the  vortex  size  and  the  boundary-layer  depth. 
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Appendix 

Derivation  of  Equation  (3) 

In  pressure  coordinates,  the  gradient  wind  equation  and 
hydrostatic  equation  may  be  written  as: 


(  3  z\ 


£  .  =  A  -  ^ 


(Al) 


and 


(A2) 


where  a  is  the  specific  volume,  p  is  the  pressure,  z  is  the 
height  of  a  pressure  surface  and  g  is  the  acceleration 
due  to  gravity.  Eliminating  the  geopotential  height  of 
the  pressure  surface,  gz,  gives  an  alternative  form  of  the 
thermal  wind  equation: 


1 


r 


3 


(A3) 


Since  s*  is  a  state  variable,  a  can  be  regarded  as  a 
function  of  p  and  5*.  Then  with  a  little  manipulation 
(14)  becomes  the  thermal  wind  equation: 


(A4) 


E86  invokes  one  of  the  Maxwell  relations  for  moist 
saturated  air  in  the  form 


so  that  (15)  becomes 


(A5) 


(A6) 


With  the  assumption  that  M  and  s*  surfaces  coincide, 
i.e.  M  —  M{s*),  (17)  becomes 


^  .  (A7) 

r3  V  3 P  Jr  \  3 P  ) s*  dM  V  3r  J 


Note  that  ( dT /dp) s*  is  just  the  temperature  lapse  rate 
as  a  function  of  pressure  along  a  moist  adiabat.  Now 
along  an  M  surface. 


/ dM\  (  dM\ 

but^Ur)^0' 


(A8) 


so  that  the  slope  of  an  M  surface  in  (r,  p)  space  is 


dr\  _  /  dM\  //3  M\ 

dp)M  \dp)r/\dr)' 
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Combining  (A7)  and  (A9),  the  thermal  wind  equation 
(A6)  becomes 


1  (dr  2\  _  1  ( dT\  ds* 

2  =  V3pA*dM’ 


(A10) 


which  may  be  integrated  upwards  along  the  M  (or  s*) 
surface  starting  from  the  top  of  the  boundary  layer  z  =  h 
to  an  outer  radius  rout  to  give 


M 


M 


-~AT-Toufs*,poat)f  (All) 
M  AM 


Assuming  that  roat  >>  r,  and  using  the  chain  rule, 
(All)  gives 


*  ds*  1 
—  Ub  -  Tout(s  ,  Pom)]  —  = 

ar  2  r 


dM1 


dr 


at  z  =  h, 

(A12) 

where  TVl  is  the  temperature  at  the  top  of  the  boundary 
layer  and  Tout  is  the  outflow  temperature  along  the  M 
(or  5*)  surface  at  rout.  Using  the  Exner  function,  n  = 
(p/Po)K,  instead  of  pressure,  the  gradient  wind  equation 
(Al)  takes  the  form 


M2  =  r3 


,  31n;r\  1  / 

CpTB^^r)z  +  4rf  _ 


(A13) 


In  the  expression  for  n,  k  =  R/cp,  where  R  is  the 
specific  gas  constant  and  pD  is  a  constant  pressure,  taken 
by  E86  to  be  1015  mb.  Substituting  (A13)  into  (A  1 2) 
results  in 


_  TB-Tout(s*,poat)d  \n6* 

TB  dr 

3  In  7r  13/  3  In  n  \  1  rf 2 

—  - 1 - (  r -  )  H - - — ,  at  z  —  h , 

dr  2  dr  \  dr  )  2  cpTB 

(A14) 

where  it  is  assumed  that  9e  =  9*  at  z  =  h.  This  equation 
is  integrated  with  respect  to  radius  from  r  to  some  large 
radius  r  =  r„,  where  it  is  assumed  that  In  n /n0  and  its 
radial  derivative  vanish,  nc>  being  the  value  of  n  at 
Z  —  h  and  r  —  ra.  Remembering  that  7b  is  assumed  to 
be  constant,  the  result  is 

*  *  1  r°  *  31n<9e 

-  In  0*  +  In  0:  +  —  Toat(s  ,  poat)  — -dr 

7b  Jr  dr 

1/  3  In  7r  \  1/  3  In  7r 

=  In  n „  —  In  jt  H —  I  r -  1 - 1  r - 

2  \  dr  ) g  2  \  dr 

1  f 2 

+  t  —  Y2  -  r2),  at  z  =  h.  (A15) 

^  C  p  1 B 
Emanuel  defines 

Y°at  =  In  (9*/9*  )  n  T°at  d(ln  )’  (Al  6) 

m(6'e/6'eo)  J\aeco 


which  is  an  average  outflow  temperature  weighted  with 
the  saturation  moist  entropy  of  the  outflow  angular 
momentum  surfaces.  Remember  that  9*  along  angular 
momentum  surfaces  is  taken  equal  to  the  equivalent 
potential  temperature,  9e,  where  the  surfaces  meet  the  top 
of  the  boundary  layer.  Then  (A15)  gives  Equation  (2). 

It  is  at  this  point  that  boundary-layer  considerations 
are  invoked.  Assuming  a  slab  boundary-layer  model  with 
uniform  depth  as  described  in  section  3.2,  E86  derives  a 
further  relationship  between  the  specific  moist  entropy  of 
the  boundary  layer,  s,  and  M  by  effectively  dividing  (10) 
by  (9).  We  recognize  here  that  the  near-surface  wind  may 
be  different  from  that  at  the  top  of  the  boundary  layer. 
Thus  following  E86,  but  allowing  for  a  reduced  surface 
wind,  we  obtain 


ds* 

AM 


z=h 


Js_ 

rM 
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(A17) 


where  xs  =  -cpCk  |VJ  (In  9e  -  ln0*)  and  rM  = 
— CoIVsIrVs  are  the  surface  fluxes  of  enthalpy  and 
momentum  expressed  by  standard  aerodynamic  formulae, 
and  |  Vs |  is  the  magnitude  of  the  near-surface  horizon¬ 
tal  velocity.  Other  quantities  are  defined  in  section  2. 
In  the  derivation  of  (A17),  it  is  assumed  that  the  spe¬ 
cific  entropy,  s,  and  the  equivalent  potential  temperature, 
0e,  are  uniform  across  the  boundary  layer  and  that  the 
air  at  the  top  of  the  sub-cloud  layer  is  saturated  so  that 
.vb  =  5*  and  9e  —  9*.  This  equation  can  then  be  blended 
with  (All)  above.  (A12)  then  gives 


In  9*  =  In  01 


at  z  —  h, 


CD 


1 


Ck  c„(Tb  —  Tout) 


V2  +  \rfV 
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where  M  has  been  expressed  in  terms  of  the  tangential 
wind  speed  V  at  z  —  h.  In  Region  II,  rf  <<  V  so  that 
the  second  term  in  parentheses  on  the  right  of  (A  18)  can 
be  neglected  compared  with  V2  and  the  equation  may  be 
written  as 


pV2  =  ^ cp(Tb  -  Tout) (In 0*  -  ln0e*),  at  z  =  h, 
to 

(A  19) 

where  p.  —  V S/V .  (A19)  is  (3)  in  section  3.1  and  is  a 
cornerstone  of  the  current  EPI-theory  (Emanuel,  1995; 
Bister  and  Emanuel.  1998). 

A  further  important  relationship  in  Emanuel’s  theory 
is  that  between  9*  and  the  pressure  and  humidity  at  the 
top  of  the  surface  layer,  which  may  be  written 


I  0e  ,  ^ 

In  —  =  —  In  — 

0*  It, a 


LqfRHA 
RY  ) 


+  -^r(RH-RHX  at  z  =  h, 

TV  i  s 


(A20) 


where  L  is  the  latent  heat  of  vaporization,  q  is  the  water 
vapour  mixing  ratio,  RH  is  the  relative  humidity,  and 
T  is  the  absolute  temperature.  The  subscript  ‘a"  denotes 


Copyright  ©  2008  Royal  Meteorological  Society 


Q.  J.  R.  Meteorol.  Soc.  134:  551-561  (2008) 
DOI:  10. 1002/qj 


HURRICANE  MODEL  AND  POTENTIAL  INTENSITY 


561 


the  ambient  value;  as  above,  a  subscript  ‘s’  denotes  a 
value  at  the  top  of  the  surface  layer  and  a  superscript  * 
denotes  a  saturation  value.  This  equation  is  the  same  as 
Equation  (25)  in  E86  if  one  assumes  that  the  reference 
pressure  in  the  definition  of  the  Exner  function  is  /;a 
rather  than  1000  mb  as  is  usual. 
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